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Synopsis 

C++ is an efficient object-oriented language of rapidly growing popularity. It can be of real value in a 
wide range of disciplines, including numerical computing, where it seems to offer important advantages over 
most competing languages. 

Object-oriented languages 

« y Vhat e * act , ly 1S an ol) j ec t-oriented language? The most important defining characteristic is support for 
polymorphic data types.” Procedural languages, like Fortran and C, contain built-in types such as integers 
reals, characters and so on. The integer type, for example, Consists of the requisite bits of data, a set of 
associated operations, +,*,/, ..., and coercions to and from the other built-in types. One can build data 
structures of arbitrary complexity in Fortran, but these are not “first class” types, like integers. 

For example, one can form a “sparse.matrix” from arrays of integers indexing into arrays of reals But 
Fortran 77 does not let one declare several of these as 

sparse .matrix A,B,C 

and then perform operations such as: 

A = B + C 

Languages like Clu and Ada, supporting “abstract data types,” let one do precisely this. One can for 
example, in Ada define a “set.oLwords” abstract data type. This would be a user defined type which might 
be useful in comparing documents. Once the type is defined, one can then declare several such sets 

set_of_words A,B,C 

One can also operate on them just as with the built in types 
A := B .+. C 

where .4*. might be a user-defined union operation. 

OO languages push this concept further, allowing one to define a “set.of.ctype T>”, where T can be 
any type in the language. This new type, a “set_of.<type T>”, is “first class” in OO languages, one can use 
variables of that type exactly like those of the built-in types. To make this clear, types are called “classes” 
in the OO world, while values of those types (classes) are called “objects,” though whehter these new terms 
do more to clarify or obfuscate is not clear. 

To see how OO ideas might be used in numerical computing, it might, for example, be useful to define 
a class “mesh-cell” which would be the basic unit of an unstructured mesh. Mesh cells come in a number 

of varieties, which can be thought of as subtypes (subclasses) of the type (class) “mesh.cell” as shown in 
Figure . 

All mesh cells share certain properties, volume, temperature, pressure etc. declared as part of class 
“mesh.cell.” Cubes and tetrahedrons share these properties, but have their own unique properties as well 
They have different numbers of faces and vertices for example. 

The ability to allow useful computing on a set of related but not identical user-defined types is the defining 
characteristic of 00 languages. In the above case, one can make an array of “mesh.c.ells” , consisting of prisms 
tetrahedrons, and cubes. One can access the volume of any element of this array, since all “mesh cells”’ 

have volume. To access specialized properties, one may have to select on the particular subclass of each 
“mesh-cell” . 
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mesh_cell 



prism tetrahedron cube 


Figure 1: Mesh cell type hierarchy 


CH — b in numerical computing 

How useful will C++ be in numerical computing? C++ contains most of the useful new features in Ada 
or Fortran 90, and is easily extensible in a number of ways. People around the world are rapidly developing 
class libraries for finite element analysis, for sparse matrix arithmetic, and so on. C++ together with a new 
class library is essentially a new application-specific, language, and one that may have a powerful impact on 
a particular subdiscipline. 

To see how this could have an impact, one need only realize that there are, for example, at least a dozen 
different unstructured grid codes here at Langley, with relatively little code shared between them. Given the 
appropriate class library supporting unstructured grids, one should be able to prototype new unstructured 
grid algorithms much faster, by borrowing large chunks of previously written code. This is the promise of 
00 computing in C++. Efficient execution, compatibility with previously written C and Fortran, and the 
00 approach are the major advantages to C++. 

C++ also has its problems. One is that its syntax and semantics, inherited from 0, are needlessly 
complex, significantly steeping the learning curve for new programmers. Another problem is that, like Ada 
and Fortran 90, C++ is a large language, full of complexities most programmers will never master. Only 
experts will master the full language, with most programmers limping along on their own particular subset. 

These problems are real, but clearly not fatal, given the exponential growth of C++. From one perspec- 
tive, C++ is essentially a halfway point between traditional procedural languages, like C and Fortran, and 
‘rapid prototyping” languages like Smalltalk. Over the longer term, as computer power increases and our 
algorithms become more complex, one expects research numerical computing, like that done at Langley, to 
shift in the rapid prototyping” direction. Use of C++ is an important step in that direction. 
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An object oriented language is one allowing users to 
create a set of related types and then intermix and ma- 
nipulate values of these related types. 


• Such types are called classes. 

• Values of such types are called objects. 

Intermixing related user-defined types is called poly- 
morphism. Polymorphism, and the things that go with 
it, encapsulation , inheritance , and dynamic binding give 
00 languages their increased semantic power. 
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Fortran is not 00. Integer variables, for example, have 
powerful properties the user cannot duplicate in new 
types. There is no way to define a sparse-matrix type, 
then do: 


sparse_matrix A,B,C 


cn 


A 


B + C 



•jo^iodo uoiun panyop i9sn v si 9i9ipw 


D •+• 9 =: V 

0‘g‘V spjo^i” jo~q.es 

: 0(Kq a spiOM~jo~q0s„ ie jo ^iijeiirosiisds,, 
v oqi^io treo 9 iio ‘ispy uj *s9dXq M9u 9uij9p oq 
sj 9 sn Avopig u s 9 di{q leq^gp qoigiqsqig,, Sniqioddns S9S'enSrn3q 


525 



526 


Typical class lattice: 


window 



bordered window 




bordered slider window 
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What is C++ ? 


• A statically- checked object oriented language down- 
ward compatible with C 

• Useful on any program with complex data struc- 
tures (though originally intended for systems pro- 
gramming) 

• Potentially useful in scientific programming, since it 
is efficient and supports complex data structures well 


Templates - data abstractions with types as parameters 

template<class T> 
class stack { 

T* v; 

T* p; 

int sz; 
public : 

stack (int s){ y = p = new T[sz = s] ; > 
"'stack () { delete [] v; } 

void push(T a){ *p++=a; } 

T pop(){ return * — p; } 

int size () const { return p-v; } 
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Operator Overloading - the ability to define new uses 
for the built-in operator symbols (+ - / ! = . . .) 

Example: polygon union 

P4 = PI + P2 + P3 

P4 - union(union(Pl ,P2) ,P3) 

P4 = PI .union(P2 .union(P3) ) 





Using inheritance in numerical codes 



cell : 

- volume 

- flow variables 

cube : 

- 8 vertices 

- 6 faces 

prism: 

- 6 vertices 

- 5 faces 
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Efficiency 


• C++ provides much of the semantic power of dynami- 
cally typed languages, like Smalltalk, but is much more 
efficient. 

• C++ code runs about as fast as C or Fortran, if one 
avoids polymorphism and abstractions. 

• The overhead in using large objects (e.g. a sparse matrix 
object) is minor. 

• Fine grained objects, such as complex variables, are inef- 
ficient and should usually be avoided. 
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Observations 


• Efficiency should be thought of in terms of the entire 
programming, debugging, and execution cycle. 

• If a language would make programming substantially 
easier, there could be such a gain from algorithmic 
improvements, that significant run-time inefficiency 
could be tolerated. 
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Conclusions 


• C++ is an effective 00 language, and has become 
the defacto standard. 

• C++ supports the data structures needed for com- 
plex numerical algorithms very well. 

• It is efficient, and can be readily intermixed with C, 
Fortran and perhaps HPF. 

• Developing numerical algorithms in C++ should in- 
crease opportunities for code reuse and for sharing 
code between programmers. 



